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ABSTRACT 

The interaction of high-velocity neutron stars with the interstellar medium produces bow shock neb- 
ulae, where the relativistic neutron star wind is confined by ram pressure. We present multi-wavelength 
observations of the Guitar Nebula, including narrow-band Ha imaging with HST/WFPC2, which re- 
solves the head of the bow shock. The HST observations are used to fit for the inclination of the pulsar 
velocity vector to the line of sight, and to determine the combination of spindown energy loss, velocity, 
and ambient density that sets the scale of the bow shock. We find that the velocity vector is most 
likely in the plane of the sky. We use the Guitar Nebula and other observed neutron star bow shocks to 
test scaling laws for their size and Ha emission, discuss their prevalence, and present criteria for their 
detectability in targeted searches. The set of Ha bow shocks shows remarkable consistency, in spite of 
the expected variation in ambient densities and orientations. Together, they support the assumption 
that a pulsar's spindown energy losses are carried away by a relativistic wind that is indistinguishable 
from being isotropic. Comparison of Ha bow shocks with X-ray and nonthermal, radio-synchrotron bow 
shocks produced by neutron stars indicates that the overall shape and scaling is consistent with the same 
physics. It also appears that nonthermal radio emission and Ha emission are mutually exclusive in the 
known objects and perhaps in all objects. 

Subject headings: starsmeutron — pulsars:individual (PSR B2224-I-65) — shock waves — ISM:general 



1. INTRODUCTION 

Bow shocks are observed on a wide variety of astrophys- 
ical scales, ranging from planetary magnetospheres (Spre- 
iter & Stahara 1995) and OB-runaway stars (van Buren & 
McCray 1988) to merging galaxy clusters (Markevitch et 
al. 2002). They have been invoked in models for beamed 
gamma-ray bursts (Wang & Loeb 2001), protostellar out- 
fiows (Ostriker et al. 2001), ultra-compact H II regions 
(van Buren & Mac Low 1992) and the interaction of the 
solar wind with the interstellar medium (Baranov, Kras- 
nobaev & Kulikovskii 1971). Some of the most spectacular 
bow shock nebulae are those associated with neutron stars. 
The spindown of neutron stars (NS) transfers rotational ki- 
netic energy into the interstellar medium, and bow shock 
nebulae provide a way to probe these energetic environ- 
ments, revealing the properties of both the NS relativistic 
winds and the ambient interstellar medium (ISM). 

The Guitar Nebula is one such visually striking bow 
shock nebula, produced by the radio pulsar B2224-I-65. 
The nebula was discovered in Ha observations using the 
5-m Hale Telescope at Palomar Observatory (Cordes, Ro- 
mani & Lundgren 1993). Until recently, there were only 
two other NS with bow shocks observed in Ha: B1957-)-20 
(Kulkarni & Hester 1988) and B0437-4715 (Bell et al. 
1995). However, Ha bow shocks have recently been dis- 
covered around a radio quiet NS, RX J1856.5— 3754 (van 
Kerkwijk & Kulkarni 2001), as well as another ordinary ra- 
dio pulsar, B0740-28 (Jones, Stappers & Gaensler 2002). 
More such nebulae will no doubt be discovered as targeted 
searches become more sensitive and cover more objects. 
Additionally, there are some bow shock nebulae which are 
detected at radio/X-ray wavelengths, but not in Ha, as- 
sociated with known radio pulsars as well as radio-quiet 
NS: for example, the "Duck" (Frail & Kulkarni 1991) is 
associated with the pulsar B 1757— 24 and the supernova 



remnant (SNR) G5.4— 1.2, while the cometary nebula in 
the SNR IC 443 is produced by a radio quiet NS (Olbert 
et al. 2001). 

In this paper, we derive scaling relationships for bow 
shock parameters based on the underlying physical pro- 
cesses, and discuss the requirements for such nebulae to 
be detectable. A detailed analysis of the Guitar Nebula is 
presented based on multi-wavelength observations, includ- 
ing high resolution Hubble Space Telescope (HST) obser- 
vations which resolve the fine structure of the bow shock. 
Through model fitting, constraints are obtained on the 
density of the ambient ISM and the inclination of the NS 
velocity vector to the line of sight (LOS). The published 
data on NS bow shocks arc consolidated with the parame- 
ters for the Guitar Nebula in order to test the derived scal- 
ing relationships for Ha bow shocks. We verify that the 
stand-off angle for plerionic bow shock nebulae detected at 
radio and X-ray wavelengths also scales in the same man- 
ner. In testing these scaling laws, it is apparent that the 
current state of available information is inadequate, espe- 
cially for radio-detected nebulae, leading to constraints on 
the ISM that are not very strong. We identify the sources 
of uncertainty, and suggest future observations that will 
be able to refine current constraints as well as confirm or 
refute physical scenarios relating to the formation of Ha 
and plerionic bow shock nebulae. 

The paper is organized as follows: the expected scaling 
relationships for bow shock parameters are derived in §2, 
along with an analysis of the detectability criteria for bow 
shock nebulae. Multi-wavelength observations of the Gui- 
tar Nebula are described in §3. Since the bright head of 
the Guitar Nebula is not resolved by ground-based optical 
images, high resolution HST observations are required to 
obtain bow shock parameters; these observations, and the 
associated model fitting, are described in §4. In §5, the 
derived scaling laws are tested for known NS bow shocks. 
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with consohdated parameters for both Ha and radio/X- 
ray bow shocks. Finally, in §6, wc summarize our results, 
discuss their implications and identify future lines of in- 
quiry. 

2. SCALING LAWS FOR BOW SHOCK NEBULAE 

Bow shock nebulae are produced when NS move super- 

sonically through the ambient ISM, creating a shocked 
layer where ram pressure balance is established between 
the relativistic NS wind and the medium. As such, bow 
shock characteristics such as the stand-off angle between 
the NS and the tip of the bow shock should scale with 
the spindown energy loss rate of the NS and the density 
of the medium. Pulsars with detected Ha bow shocks 
(J0437-4715, B0740-28, B1957-h20 and B2224-F65) have 
a wide range of properties, with E ranging over two or- 
ders of magnitude, and velocities and distances spanning 
an order of magnitude, as detailed in §5. This leads to the 
question of how the bow shock scales with changes in the 
properties of the NS and the ambient medium. 

2.1. Scaling of the Stand- Off Angle 

Wilkin (1996) provides an exact analytic solution for a 
stellar wind bow shock in the thin-shell limit, a problem 
first solved numerically by Baranov, Krasnobaev & Ku- 
likovskii (1971). In this analytic treatment, ram pressure 
balance between the stellar wind and the ambient medium 
determines the stand-off radius of the bow shock (i?o). For 
a star moving with velocity v^, through a uniform medium 
of density pa, with a stellar wind mass loss rate rhw and 
wind velocity v^, the condition for ram-pressure balance 
is 

leading to an expression for the stand-off radius, 

/ • \ 
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For a relativistic NS wind, the momentum outflow per 
unit area rtiwVw can be recast as E/c under the assump- 
tion that the spindown energy loss is carried away by the 
relativistic wind. Then the stand-off radius becomes 

^0 = (tz^I • (3) 
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The substitution of a relativistic NS wind in place of an 
isotropic stellar wind may be problematic, since the NS 
wind is at best only quasi-isotropic. Some authors have 
cited the possibility of an aligninc;nt of the NS spin axis and 
proper motion (e.g. Spruit & Phinney 1998), and two pul- 
sars (the Crab and Vela) apparently conform to this pic- 
ture (Lai, Chernoff & Cordes 2001), although it might not 
be a general phenomenon (Deshpande, Ramachandran & 
Radhakrishnan 1999). In case of alignment, the rotation- 
averaged relativistic wind, while not isotropic, would be 
axisymmetric about the velocity vector, a situation which 
adds only one extra scaling parameter to the bow shock 
description (Wilkin 2000). If there is no such alignment, 
so that the rotation vector is skewed with respect to both 
the velocity and the magnetic field, the behavior of the 
relativistic wind depends on the relative orientations of 



these vectors and the opening angle of the wind: for a 
range of orientations, rotational averaging leads to a rjuasi- 
isotropic situation relative to the velocity vector. Specifi- 
cally, even for B2224-I-65, a relatively slow-spinning pulsar 
with a small stand-off radius Rq, the light cylinder radius 
(at which point the rotation-averaged NS wind has spread 
out over a large solid angle) is ~ a few xlO~^Ro. In the 
absence of additional information, we adopt an isotropic 
description for the relativistic wind. 

The observed stand-off angle Oa is simply the projected 
stand-off radius at the NS distance: 

60 = Rosmi/D, (4) 

where i is the inclination angle to the LOS, so that i = w/2 

corresponds to a velocity vector in the plane of the sky, and 
i — corresponds to motion towards the observer. The 
measurable transverse velocity of the NS, v±, is related to 
the actual velocity by the same projection factor: 

v± = fjtD = Vi, sini, (5) 

where /z is the measured angular proper motion. The am- 
bient density of the ISM, pa, is the product of the number 
density of the medium (nu), the mass of the H atom (toh) 
and the equivalent molecular weight of the hydrogen and 
helium mixture (7H = 1.37 for cosmic abundances, with 
27% He by mass). Defining riA, the density in atomic mass 
units per unit volume. 

Pa = riHlHmH = TlAmH- (6) 

Making the appropriate substitutions in Eqn. 3, and fac- 
toring in terms of constants, unknown quantities and ob- 
servables, we obtain: 
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the last factor comprises the quantities E, D and /i, which 
are (at least in principle) directly measurable. Eqn. 7 can 
be evaluated numerically, for ^0 in milliarcsecond, £^33 in 
10^3 erg s~i, D in kpc, /xioo in 100 milliarcsecond yr"^ 
and UA in cm""^: 
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2.2. Scaling of the Ha Flux 



Bow shocks emit Ha photons when the partially neutral 
ISM encounters the shock front, traveling at the NS veloc- 
ity v^,. The neutral atoms are not immediately affected 
by the shock transition, but undergo collisional excitation 
and ionization in the hot post-shock flow due to encounters 
with shocked electrons and protons (Chevalier & Raymond 
1978). Ha emission is proportional to the probability of 
an excitation before ionization, which is given by the ratio 
of the excitation rate ^ex to the ionization rate (Ray- 
mond 1991). The observed rate of Ha emission (Fa) from 
a nebula which subtends a solid angle Af2 is given by: 

Fry = — nuiv* , (9) 

9i 47r5fHa 



where gjja is the extinction correction for Ha (w 3 for the 
Guitar Nebula), and Raymond (1991) gives Qex/qi ~ 0.2 
Ha photons per neutral atom. 
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The solid angle subtended by the bow shock nebula 
Af2 ^ (rjB.o/D)'^, where 77 is a geometric factor which de- 
pends on the NS velocity. Cordes, Romani & Lundgren 
(1993) use T] oc v^,, which is physically plausible, since a 
faster-moving NS sweeps up a larger volume of ISM within 
the shock front during the ionization timescale, so that 
AO oc w^. To retain generality, we use rj ocvt, which im- 
plies AJ7 cx v1^{Rq/DY. In addition, defining the neutral 
fraction X = nni/n^ leads to riHi = Xn^i oc Xpa, and 
substituting for Af2 and riHi in Eqn. 9, 

F^o^Xnnv^inRo/Df. (10) 

The Ha photons are emitted isotropically by coUisionally 
excited neutral atoms, so that geometric projection of the 
stand-off radius and the nebula itself does not affect the 
total observed flux Fa, but projection effects have to be 
accounted for to express the actual NS velocity t;* in terms 
of the observables /i and D. Thus the Ha flux can be 
expressed in terms of observable quantities {9o,l^,D) and 
unknowns (sini, rin, -^): 

/ ^ X 2/3+1 

Fa(xXnKelvy+^ (xXnaeU^) . (11) 

\smi J 

Alternatively, using the definition of Rq from Eqn. 3, 
a scaling relationship can be derived for the observed Ha 
flux in terms of measurable quantities and the unknown 
neutral fraction: 

,,2/3-1 
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which reduces to the expression derived by Cordes, Ro- 
mani & Lundgren (1993) for = 1. 

2.3. Detectability of Ha Bow Shock Nebulae 

For bow shock nebulae, the observed Ha emission comes 
from the outer shocked layer, where the neutral medium is 
swept up and coUisionally excited. There is also an inner 
shocked region which encloses the relativistic NS wind and 
may produce synchrotron radiation detectable in X-ray or 
radio observations. Ha emission reqiiires the existence of 
a significant neutral component in the medium, which is 
true for only - 20% of the ISM (McKee & Ostriker 1977). 
Thus the non-detection of a nebula in Ha may reflect pre- 
ionization of the ISM by supernovae, or thermal radiation 
from the NS or shocked gas, and places an upper limit on 
the neutral fraction present in the ambient medium. This 
is especially applicable for nebiflae detected at radio/X-ray 
wavelengths, but not seen in Ha. 

Detection criteria can be formulated in order to predict 
whether a pulsar should produce a bow shock nebula visi- 
ble in Ha. Considering the average flux per detector pixel, 
the ratio of the Ha flux to the angular size of the nebula. 
Fa/ Aft, must exceed a (detector dependent) threshold: 

AT " = AO /A 'S'thresh- (13) 

-'Vpixel AS2/AcCD 

From Eqn. 9, this implies that Xrinv^,/ gna must exceed a 
limiting value, a result independent of the spindown flux 
E/D"^. However, a criterion based on the average flux 
does not account for brightness variations within the neb- 
ula, which can be substantial (as seen, for example, in the 
Guitar Nebula, discussed in §3). The optimal detection 



method is to block average all pixels which contain nebu- 
lar emission, in which case (assuming uncorrclatcd random 
noise in each pixel) the average flux per detector pixel only 
needs to exceed ^thresh/ v^-^pixei for a detection: 



Fr, 



Nr. 



> 7'"°^'^ (14) 

Vpixel Y iVpixel 

It is worth noting, however, that the optimal averaging 
size may not be known in advance, and a computationally 
intensive search may be required. For optimal detection, 
again using Eqn. 9 and substituting for the angular size of 
the nebula, 
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where k encapsulates various constants and numerical fac- 
tors. As expected. Ha bow shock nebulae are more likely 
to be detected for pulsars with larger spindown flux and 
higher velocities, located in regions with higher neutral 
fractions and lower extinction. 

The detection criterion formulated in Eqn. 15 can be 
recast in terms of the minimum NS velocity required to 
produce a detectable bowshock, given the spindown flux 
and the ambient ISM. For concreteness, we assume a tele- 
scope sensitivity comparable to the 5-m Hale Telescope at 
Palomar, as used to detect the Guitar Nebula (§3), and 
/3 = 1 as before. Then, for D in kpc, E in units of lO'^'^ 
erg s~^, and riH in cm~^, 

-1 gaa / ^^33 



14 km s 



Xn 



1/2 



(16) 



The relationship derived for the minimum required ve- 
locity is plotted in Fig. 1, along with transverse velocities 
for all pulsars with measured proper motions. All known 
bow shock nebulae are marked on the plot: their locations 
confirm the applicability of the detection criterion devel- 
oped above. As mentioned above, the X-ray and radio- 
detected bow shock nebulae which are not detected in Ha 
provide upper limits on the neutral fraction present in the 
ambient medium. The plot also identifies future candi- 
dates for deep searches for bow shock nebulae, although 
the unc-ertaiiities in the distance and velocity estimates 
for most pulsars and the large range of plausible values 

1/2 

for Xn^ /gua are likely to limit its usefulness as a strong 
predicting tool. 

3. MULTI-WAVELENGTH OBSERVATIONS OF THE GUITAR 

NEBULA 

The visually spectacular Guitar Nebula is a bow shock 
nebula produced by the otherwise unexceptional pulsar 
B2224-F65 (P = 0.68 s, P = 9.7 x IQ-^^ s s^^). The 
spindown luminosity of the pulsar {E = IQ-^-^ "® erg s~^, 
assuming a moment of inertia of 10^^ gm cm^) is sig- 
nificantly lower than the E of other pulsars with known 
Ha bow shocks. However, B2224-I-65 has an extremely 
large space velocity, which evidently compensates for the 
low E in producing a visible bow shock nebula. Its 
proper motion = 182 ± 3 mas yr^^, at a position an- 
gle 52.1° ± 0.9° (Harrison, Lyne & Anderson 1993). The 
dispersion measure (DM = 36.16 ± 0.05 pc cm~'^) of the 
pulsar, combined with a model for the Galactic electron 
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Fig. 1. — Detectability of bow shock nebulae. Transverse veloci- 
ties are plotted against the spindown flux E/D'^ for all pulsars with 
published proper motions. The velocity and distance estimates are 
significantly better for pulsars with parallax distances, which are 
plotted with larger open circles. Lines are overplotted for the de- 

1/2 

tcction criterion in Eqn. 16, for Xn-^ /9a = 0.01 (dashed) and 0.1 
(dotted). All known H« and radio / X-ray bow shock nebulae are 
identified on the plot. The Crab (B0531+21) and Vela (B0833-45) 
pulsars are not detected as bow shock nebulae because of the negli- 
gible neutral fraction in their local ISM. 



density distribution (Taylor & Cordes 1993) yields a dis- 
tance £) « 2.0 ± 0.5 kpc, at which the transverse velocity 
v± = Dfj, w (1.7 ± 0.4) X lO'' km s"""^, the largest known 
velocity among radio pulsars. 

It is evident from the discovery images that there may 
be a significant amount of ionized material in the vicinity 
of the pulsar, biasing the DM upward, and hence the dis- 
tance may be overestimated. However, even at a distance 
of 1 kpc, the pulsar velocity (> 850 km s^^) is significantly 
higher than the mean population velocity of ^ 450 km s^^ 
estimated by Lyne & Lorimer (1994). Recent studies (Ar- 
zoumanian, Chernoff & Cordes 2001; Cordes & Chernoff 
1998) prefer a bimodal velocity distribution, with a signif- 
icant portion (^ 15%) of pulsars having a space velocity 
> 1000 km s"^ B2224+65 is probably at the high end of 
this distribution, and the bow shock nebula is visible in 
spite of the small E due to the large space velocity of the 
pulsar. 

3.1. Ground-based Optical Observations 

In an ongoing program at Palomar Observatory, the 
Guitar Nebula has been imaged regularly at the 5-m Hale 
Telescope using the Carnegie Observatories Spectroscopic 
Multislit and Imaging Camera (COSMIC; Kells et al. 
1998), mounted at the f/3.5 prime focus, in reimaging 
mode. A narrow- band Ha filter (20 A at 6564 A) was 
used. The typical seeing was 0.9"-1.4", significantly 
larger than the CCD pixel scale (0.4"). The CCD data 
were processed in standard ways with IRAF tasks. Data 
obtained in 1995 were combined (9000 s total integration 
time) to produce the representative image presented in 
Fig. 2. 

The Palomar image of the Guitar Nebula shows a re- 



markable structure consisting of a faint, limb-brightened 
"guitar body" which is approximately symmetric, with an 
elongated neck and a bright head, the tip of which co- 
incides with the location of the pulsar. The axis of the 
nebula matches the position angle of the pulsar proper 
motion (52.1° ± 0.9°), as expected. Differential Galactic 
rotation causes only a small effect (< 3% of the transverse 
velocity) and is therefore neglected. 

The body of the guitar appears to be composed of two 
closed-off bubbles. This could be explained by the pro- 
jection of an open structure onto the plane of the sky, 
or by episodic variation in E, which is not inconceivable 
given that B2224-I-65 has been known to glitch. The com- 
plex shape may alternatively be explained by the pene- 
tration of neutral ISM atoms into the space occupied by 
the pulsar wind, as suggested by Bucciantini & Bandiera 
(2001). However, the simplest explanation is the existence 
of density variations in the ambient medium. The observed 
closed bubble structures can be produced by density gra- 
dients oriented roughly parallel to the NS transverse ve- 
locity. The nebula brightens where it narrows between the 
two bubbles, as expected in a higher density region. The 
extremely bright head is then created by the NS entering 
a region with much higher (~ 10 x) density, or a signifi- 
cantly higher neutral fraction. The distinct variations in 
brightness along the nebula, the observed departures from 
symmetry on small scales and the existence of large-scale 
filamentary structures (visible at the bottom of Fig. 2) 
can all be explained within this scenario, although other 
explanations cannot be ruled out. 

3.2. Radio Interferometric Observations 

The Guitar Nebula was observed at radio wavelengths 
with the NRAO Very Large Array (VLA). In order to get 
comparable resolutions over a range of frequencies, differ- 
ent VLA configurations were used: B-array for 1.4 GHz 
and 4.8 GHz, and C-array for 8.4 GHz observations. 
About 3 hours of integration time were obtained at each 
frequency, though the 8.4 GHz observation was plagued by 
poor weather conditions and did not yield as much usable 
data as the other frequencies. The observation parameters 
are summarized in Table 1. 

The data were deconvolved using the CLEAN algorithm 
(Hogbom 1974), as implemented in AIPS, the Astronomi- 
cal Image Processing System, in order to remove the side- 
lobes of strong sources within the primary beam. Wide- 
field images were produced at each frequency with map 
noise ranging from 23 /uJy to 28 yuJy, as detailed in the 
last column of Table 1. 

The radio pulsar B2224-I-65 is detected at 1.4 GHz as 
a point-source with a flux ~ 3.2 mjy. It is luidctected 
at the higher frequencies. No extended features related 
to the Guitar Nebula or the optical filamentary struc- 
tures are observed. We derive an upper limit (^ 5(t) 
of 0.1 mJy beam^^ on any synchrotron emission, either 
from the relativistic pulsar wind, or from the shock accel- 
erated particles in the nebula, at wavelengths between 1.4 
and 8.4 GHz with 1-4" beam sizes. The non-detections 
thus rule out any plerionic component to the pulsar or 
any non-thermal emission from the Guitar Nebula over 
the observed range of wavelengths and resolutions. The 
data were also smoothed by tapering in the U — V do- 
main: although the possibility of faint extended emission 
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Fig. 2.— The Guitar Nebula in Ha (20 A filter at 6564 A), imaged with the 5-m Hale Telescope at Palomar Observatory in mid-1995. 
North is up and east is to the left, as indicated, and the L has arms 10" in size. 



Table 1 

Multi-array Radio Observations at the VLA 



t'obs 


VLA 


Beam 


Tint 


Az/ 


a (RMS) 


(GHz) 


Array 
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(hr) 


(MHz) 


(^Jy/bcam) 


1.4 


B 


3.8 


2.9 


50 


27 


4.9 


B 


1.2 


2.9 


50 


23 


8.4 


C 


2.2 


2.4 


100 


28 



on much larger scales cannot be excluded, no significant 
flux was detected from the Guitar Nebula after smoothing 
to 2-8" beam sizes, comparable to the size of the bright 
part of the nebula in Ha. 

The derived upper limit on plerionic radio emission is 
more stringent than those derived by Gaensler et al. (2000) 
for other pulsars, and the non-detection is consistent with 
their conclusion that older pulsars may become less effi- 
cient at producing synchrotron emission with increasing 
age. 

3.3. X-ray and Infra-Red Observations 

X-ray observations of the Guitar Nebula using the 
ROSAT High Resolution Imager have been described by 
Romani, Cordes & Yadigaroglu (1997). Weak soft X-ray 
emission is detected using an aperture mask matched to 
the optical nebula, but the detection is significant at only 
^ 4(T. Follow-up observations using Chandra have recently 
been obtained, and will be described elsewhere. 

Narrow-band infrared observations with PFIRCAM at 
the Hale 5-m telescope at Palomar Observatory resulted in 
non-detections, which complement the radio observations 
in limiting any thermal or non-thermal emission from the 
Guitar Nebula. Upper limits on the infrared flux will be 
published in future work. 



4. HUBBLE SPACE TELESCOPE OBSERVATIONS AND 
MODELING 

The head of the Guitar Nebula is not well resolved 
by ground-based observations, and the expected scale of 
the shocked region where the NS wind meets the ambi- 
ent medium (^^ 0.05") cannot be probed from the ground. 
HST observations are required, followed by model fitting 
to obtain the bow shock parameters. 

4.1. HST Observations 

High-resolution HST observations were obtained in 1994 
December, using the Wide Field and Planetary Camera 
2 (WFPC2; Holtzman et al. 1995). The bright head of 
the nebula was centered on the Planetary Camera (PC), 
which provides 0.0455" resolution. The observations in- 
clude 6 dithered exposures of 1200 seconds each, using fil- 
ter F656N (22 A at 6564 A), comparable to the filter used 
at Palomar. The images were combined using Variable- 
Pixel Linear Reconstruction (or 'Drizzling'), as described 
by Fruchter & Hook (2002). This includes cosmic ray re- 
jection and partially compensates for the undersampling 
of the HST point-spread function by the 0.0996" pixels of 
the WF chips on WFPC2. 

A section of the PC image is presented in Fig. 3: for 
reference, the three stars visible in the frame are also seen 
in the Palomar image (Fig. 2). In spite of the limited 
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Fig. 3.— The head of the Guitar Nebula, imaged with the HST Planetary Camera (7200 seeonds, 22 A filter at 6564 A) in late f994. 
The left panel shows a seetion of the image with a hyperbolie-sine streteh applied to the grayseale. The lower right panel marks the region 
eovered by the image section, superimposed on the near-contemporaneous Palomar image (mid 1995). The upper right panel shows the HST 
image convolved with the typical seeing disk (~ 1" FWHM) at Palomar, which makes the bright head and the neck of the Guitar visible. For 
reference, the three bright stars in the HST images are also visible on the Palomar image (Fig. 2) . (Image degraded to fit astro-ph size limits) 



signal-to-noise ratio (S/N) of the image, it is apparent that 
the bright head of the Guitar Nebula, unresolved in the 

Palomar images, is resolved by the PC. The bow shock is 
limb-brightened, as observed in the extended Guitar body 
observed at Palomar. The bright head and the neck of 
the Guitar become visible when the PC image is smoothed 
down to the typical seeing at Palomar, as shown in the up- 
per right panel of Fig. 3. The extended body of the Guitar 
and the nearby filamentary structures are also detected in 
the WF images after smoothing (not shown here). 

Along with the narrow-band images, broad-band images 
were also obtained with WFPC2, with 4 dithered expo- 
sures of 180 seconds each using filter F675W (867 A at 
6717 A). The combined image (720 s) shows no significant 
emission from the nebula. As expected, no optical coun- 
terpart is detected for the radio pulsar, down to a limiting 
magnitude of i? = 25.8 (Scr) on the broad-band PC image. 

4.2. Evolution of the Guitar Nebula 

As discussed above, the simplest explanation for the 
characteristic shape of the Guitar Nebula is the propa- 
gation of the NS through a non-uniform ambient medium, 
in which case the shape of the nebula provides information 
about the density of the medium encountered in the past. 

Fig. 4 shows the symmetric and anti-symmetric compo- 
nents of the head of the Guitar Nebula, which were ob- 
tained by mirroring a smoothed version of the PC image 
about the axis of the nebula and taking the sum and the 
difference of the two. The axis was found by using the 
pulsar proper motion as a starting point, and then fit- 
ting for the bow shock position angle, as described below. 
The symmetric component of the image shows significant 
changes in the opening angle of the nebula at 2" scales. 



which are also visible, though less prominent, in the un- 
smoothed HST image (Fig. 3). It appears that the opening 
angle has narrowed recently (< 10 yr, between 1 and 2" 
downstream from the tip), probably as the NS encountered 
a gradient of increasing density. 




Fig. 4. — Symmetric and anti-symmetric components of the head 
of the Guitar Nebula, obtained by mirroring the head of the nebula 
about its axis. The image has been enhanced to bring out detail by 
assigning the mean amplitude to all pixels with amplitudes within 
±1(T of the image mean. The range of the greyscale is the same for 
both panels. While the tip of the head is almost perfectly symmetric, 
the limbs of the nebula are visible in the anti-symmetric image at 
the ±2<T level. 

There is also a small but significant departure from bi- 
lateral symmetry along the head of the Guitar Nebula, as 
illustrated by the image of the anti-symmetric component. 
The tip of the head is almost perfectly symmetric, and dis- 
appears in the anti-symmetric image, but the north-west 
and south-east limbs of the nebula deviate from the axis of 
symmetry at the 2a level, consistent with a density gradi- 
ent in the ambient medium which is not aligned with the 
NS velocity. 
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In both the smoothed HST image (Fig. 3) and the im- 
age of the symmetric component, the head of the nebula 
appears to have a closed-off structure, reminiscent of the 
rounded end of the "guitar body" seen on the Palomar 
image (Fig. 2). Once again, while various other explana- 
tions (including a time- variable E) cannot be ruled out, 
the simplest explanation involves the existence of a den- 
sity gradient, which confines the expansion of the head of 
the nebula. This explains both the brightening of the head 
of the nebula and the faint tubular "neck" of the guitar. 

As part of an ongoing effort to monitor the evolution 
of the nebula, several images have been obtained at the 
5-m Hale Telescope at Palomar Observatory over the past 
7 years. Additionally, we expect to obtain another HST 
observation at higher S/N in late 2001, at which time the 
NS will have moved ~ 1" from the previous position. This 
will allow a detailed study of the dynamic evolution of the 
head of the nebula and the ISM in the neighborhood of 
the NS, or alternatively, more exotic scenarios with time- 
variable E. 

4.3. Modeling the Bow Shock 

The lack of an optical counterpart to PSR B2224+65 

makes it difficult to obtain bow shock parameters like the 
stand-off angle 6q by inspection of the image. Model fitting 
must be used in order to extract from the HST image, 
which resolves the stand-off region at the head of the bow 
shock. 

Adopting the simple momentum-conserving description 
from Wilkin (1996), the shape of the bow shock can be 
expressed in polar coordinates {R, 9) with the NS at the 
origin and traveling along the z axis: 

1/2 



sin^ 



1 



tan 6* 



(17) 



This expression describes the shape of a thin shell bow 

shock provided that the shocked gas undergoes rapid ra- 
diative cooling. If the cooling is not rapid enough, ther- 
mal pressure should result in a somewhat wider tail down- 
stream along the bow shock. The momentum-conserving 
description also assumes mixing in the shocked layer be- 
tween the ambient medium and the relativistic wind, 
whereas a laminar flow may be a more accurate description 
(Bucciantini & Bandiera 2001). The model also assumes 
an ambient medium with uniform density, which, we have 
argued, is not the case here. However, the present work is 
limited to extracting basic bow shock parameters from an 
image with poor S/N. To work around these limitations, 
the model fitting procedure was restricted to within 2.6" 
of the tip of the bright head, where the opening angle does 
not change much, the surface is relatively smooth, and the 
behavior of the tail downstream is not relevant. 

The thin shell bow shock described by Eqn. 17 can be 
scaled and projected to account for viewing geometry. For 
each set of scale and geometry parameters, the projected 
surface density is calculated by integrating the density 
through the bow shock along each line of sight. The Ha 
flux from the nebula is proportional to the projected sur- 
face density, accounting for the observed limb-brightened 
appearance. 

A bow shock nebula described by this model can be pa- 
rameterized by two angles, the inclination angle of z to the 



LOS (i = for motion towards the observer, and i = 7r/2 
for motion in the plane of the sky), and a position angle 
on the sky, measured east of north. Additional parameters 
are the thickness of the shocked layer which emits in Ha 
as a fraction of the stand-off radius, a constant numerical 
factor which relates the Ha flux to the projected surface 
density, and a dimensionless scale parameter, 

S = D^^lnl'/'{E/Ejj'/\ (18) 

which accounts for the apparent change in size with dis- 
tance D, the number density of the ambient medium ha, 
and the ratio of the actual NS luminosity to the calculated 
spindown luminosity. 

Models were generated over a grid of parameter values 
for position angle (43°-53°), scale parameter S (1.0-6.0), 
inclination to the LOS i (15°-90°) and the thickness of 
the Ha-emitting layer (5% 25% of the stand off radius). 
These models were compared to the part of the Guitar 
Nebula least affected by the non-uniformity of the ambient 
medium, within ~ 2.6" of nebula tip in Fig. 3, which covers 
~ 6, 600 pixels on the PC image. Best fit parameters were 
obtained using least-squares minimization of the difference 
between the data (Dk) and the model (M^) over the image 
{1 < k < N, where A'' is the number of pixels), 

N 

=Y,a-\Dk-aMkf, (19) 

fc=i 

where a is the numerical factor relating the Ha flux to 
the projected surface density, calculated using a matched 
filtering approach: 

= . (20) 

The model fitting procedure yields a good estimate for 
the stand-off angle 9q separating the shock front and the 
NS. For the entire range of plausible models, we obtain 
9(j = 0.06" ± 0.02". The position angle is also determined 
to be 48° ±2°, in agreement with the pulsar proper motion 
to within the errors. The measured stand-off angle is used 
to test scaling laws for bow shock nebulae in the next 
section; we note here that the measured 6*0 can be used in 
Eqn. 8 to constrain the distance to B2224-I-65: 

, sin j 



pc — 0.75 



(21) 



However, model fitting was inconclusive for the other 
parameters, especially the inclination angle i. A mini- 
mum was found in the surface at i = 90°, with reduced 
= 1.20, compared to a reduced = 1-30 for a model 
consisting of a flat background only. However, the surface 
is not simple, and shows another secondary minimum at 
i ^ 30° (reduced x^ — 1.21). The corresponding values 
of the scale factor S are 5.0 and 2.5 respectively. Phys- 
ically, i ~ 90° is preferred, since a low inclination angle 
implies an extremely high three-dimensional velocity. The 
primary reason for the poor fit is the low S/N in the im- 
age. Observations scheduled for late 2001 with the HST 
WFPC2 will have a signiflcantly longer integration time, 
and hence a better S/N, which is expected to allow a more 
discriminating model fit. 

The narrowing of the opening angle of the bow shock 
and the departure from bilateral symmetry, which are il- 
lustrated in Fig. 4, may be modeled by a density gradi- 
ent in the ambient medium which is not aligned with the 
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Table 2 

Iln BoA\' Shock Parameters 



Pulsar 


D 


M 


logio E 


Oo 




Refs. 




(kpc) 


(mas yr~^) 


(erg s~i) 


( ) 


(cm s ) 




J0437-4715 


0.139 ±0.003 


140.892 ±0.006 


34.07 


10.0 ± 1.5 


2.5 X 10-3 


1,2 


B0740-28 


1.9 ±0.5 


29.0 ±0.9 


35.16 


1.0 ±0.2 


5 X 10-'"^ 


3 


B1957+20 


1.5 ±0.4 


30.4 ±0.6 


35.20 


4.0 ±0.6 


1.09 X 10-'' 


4 


B2224+65 


2.0 ±0.5 


182 ±3 


33.08 


0.06 ±0.02 


1.1 X 10-3 


5,6 


J1856.5-3754 


0.14 ±0.04 


332.0 ±0.8 


32.9 ±0.3 


1.0 ±0.2 


2 X 10-^ 


7,8,t 



References. — (1) Bell et al. (1995); (2) van Straten et al. (2001); (3) Jones, Stappers & Gaensler 
(2002); (4) Kulkarni & Hester (1988); (5) Cordes, Romani & Lundgren (1993); (6) this work; (7) 
Kaplan, van Kerkwijk & Anderson (2002); (8) van Kerkwijk & Kulkarni (2001). (\: Not detected 
as a radio or X-ray pulsar.) 



NS velocity vector. Wilkin (2000) provides a more gen- 
eral treatment than the current one, and this may be used 
for future modeling, along with a full numerical solution. 
However, the current image S/N does not justify more 
complex models with additional free parameters. 

5. TESTS OF SCALING LAWS 

5.1. Ha Bow Shocks 

The bow shock parameters obtained for the Guitar Neb- 
ula can be used, in conjunction with parameters measured 
for other known bow shock nebulae, to test the scaling 
laws derived in §2. Besides the Guitar Nebula, bow shocks 
have been detected in Ha for the radio pulsars B1957+20 
(Kulkarni & Hester 1988), B0437-4715 (Bell et al. 1995) 
and B0740-28 (Jones, Stappers & Gaensler 2002). In each 
case, the spectra are dominated by Hydrogen Balmer lines. 
While the proper motion and E of these pulsars are well 
measured, the distances are obtained from the pulsar DM 
along with a model for the Galactic electron density (Tay- 
lor & Cordes 1993), and are therefore relatively poorly 
constrained. The exception to this is PSR B0437-4715, 
which has an extremely precise parallax measurement from 
timing observations (van Straten et al. 2001). 

The observable parameters for these pulsars are summa- 
rized in Table 2, with ad hoc 15-20% errors assigned to 0q 
measurements except for the Guitar Nebula. It is worth 
noting that E spans more than two orders of magnitude 
for these four pulsars, and the proper motions and dis- 
tance estimates span an order of magnitude as well. The 
scaling of the stand-off angle for Ha bow shocks is demon- 
strated in Fig. 5, where the observed stand-off angle is 
plotted against the other observable quantities, and lines 
of constant are overplotted for i = 90° and 45°. Due 
to projection along the LOS, 9^) < B.q/D. Thus, for each 
observed bow shock nebula, the lines of constant n\ for 
i = 90° place an upper limit on the density of the ambient 
medium, independent of the actual inclination to the LOS. 
For example, for the Guitar Nebula, the nominal ambient 
density is constrained to be ua ^ 0.015 cm-^, which cor- 
responds to riH ^ 0.011 cm-3 for a Helium fraction of 0.27 



by mass (cosmic abundance) in the ISM (i.e. 7h = 1.37). 

The scaling laws derived for Ha bow shocks associated 
with radio pulsars can also be generalized to radio-quiet 
NS, under the assumption that these NS have similar rela- 
tivistic winds. RX J 1856. 5— 3754 is the nearest known NS, 
with a published HST parallax distance of 61 ± 9 pc (Wal- 
ter 2001). This distance has been called into question by 
Kaplan, van Kerkwijk & Anderson (2002), who reanalyze 
the same data to obtain a parallax distance = 140 ±40 pc, 
which we adopt here. While the lack of radio pulsations 
can be explained by postulating a mis-directed pulsar, the 
absence of detected X-ray pulsations from thermal hot 
spots (Ransom, Gaensler & Slane 2002) is more puzzling. 
An associated cometary Ha nebula has been detected by 
van Kerkwijk & Kulkarni (2001), who describe it as being 
consistent with either a bow shock or an ionization nebula. 
Assuming that it is a bow shock (which is the simpler of 
the two possibilities, especially in the light of the revised 
distance) and using a formalism similar to this work, van 
Kerkwijk & Kulkarni (2001) estimate i; ~ 6 x lO^^ erg s-\ 
which increases to ~ 8 x 10^^ erg for the revised dis- 
tance. From Fig. 5, tt-a ^ 1 cm-^, consistent with the 
value derived by van Kerkwijk & Kulkarni (2001). 

It is worth emphasizing the diversity of the bow shock 
nebulae in Fig. 5. B1957+20 and J0437-4715 are millisec- 
ond radio pulsars while RX J 1856. 5— 3754 is a radio-quiet 
NS, and all three produce canonical bow shock nebulae. 
On the other hand, B2224±65 and B0740-28 are ordinary 
radio pulsars but the bow shock nebulae they produce do 
not have simple shapes: the former produces the Guitar 
Nebula while the latter produces a nebula shaped like a key 
hole. Faint X-ray emission has been detected by ROSAT 
for B2224+65 (Romani, Cordes & Yadigaroglu 1997) and 
B1957±20 (Fruchter et al. 1992), although in neither case 
is the emission definitively associated with the bow shock 
nebula. In spite of this diversity, the stand-off angle, at 
least, is well parameterized by a few simple physical quan- 
tities, as expected from the underlying physics. 

The total Ha flux from a bow shock nebula also scales 
with a combination of measured, observable and unknown 
parameters, as described by Eqns. 11 and 12. While E, /x. 
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Fig. 5. — Scaling of the stand-ofF angle for Ha bow shock nebulae 
with other observable parameters, as listed in Table 2. Circles mark 
radio pulsars and the cross marks RX J1856.5— 3754, a radio-quiet 
NS. Lines of constant number density n\ = 0.01, 0.1, 1.0 cm~^ are 
overplotted for inclination to the line of sight i = 90° (thinnest to 
thickest solid lines) and i = 45° (thinnest to thickest dashed lines). 

D and 0o arc accessible to observation, the density tih is 
unknown, and few a priori constraints exist on the neutral 
fraction X. Additionally, the scaling is a function of how 
the size of the bow shock depends on the velocity, param- 
eterized as Ail cx wj^, and we have argued that (3=1. 

The five bow shock nebulae with published Ha flux val- 
ues {Fa) are plotted in Fig. 6 against two different com- 
binations of observable and measured parameters, with 
(3 = 1. While a correlation is apparent in the plot of Fa 
against 6^^^D^, it is worth noting that differences in the 
density and neutral fraction in the environment of each ob- 
ject and the different extinctions along each LOS are not 
accounted for in this plot, and neither is the projection 
effect for the velocity (sin'^i), which differs for each bow 
shock. On the other hand, for /3 = (nebula size indepen- 
dent of velocity) and /3 = 0.5, the points appear scattered 
with no discernible trend, favoring /3 = 1 and Afl cx wj. 
The two combinations of observables and E^/D) 

are related to each other through Eqn. 7, with differences 
in the medium density accounting for the relative scatter 
between the two plots. 

We note that there are two independent observed quan- 
tities (^0 and Fa) which depend on the measurable values 
for E, fi and D (as derived in §2). In principle, these can 
be combined to jointly constrain the unknown values of 
the density riH, neutral fraction X and the inclination to 
the LOS. While current observations still lack the preci- 
sion to produce meaningful constraints on all the unknown 
parameters at the same time, a parallax distance to PSR 
B2224+65 or a more precise measurement of the bow shock 
parameters for J0437— 4715, for example, can be exploited 
in the future to obtain such constraints. 
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Fig. 6. — Scaling of Ha flux from bow shock nebulae with two 
different combinations of observable parameters, for E in erg s"'^, 
9o in arcsec, /i in mas yr~^ and D in kpc. Circles mark radio pulsars 
and the cross marks RX J1856.5— 3754, a radio-quiet NS. The effects 
of the different inclinations to the LOS and the differing extinctions, 
densities and neutral fractions are not accounted for in these plots, 
and axe responsible for the relative scatter between the two panels. 

5.2. Bow Shocks at Radio/X-ray Wavelengths 

As mentioned in §2.3, NS bow shocks consist of two 
shocked regions: an outer shock, where the ambient ISM 
undergoes collisional excitation and may emit in Ha, and 
an inner shocked region enclosing the relativistic NS wind, 
which may emit synchrotron radiation. No bow shocks 
have currently been detected in both Ha and synchrotron 
radiation at radio wavelengths, although a faint extended 
X-ray counterpart may exist for the Guitar Nebula. It 
is possible that synchrotron-bright nebulae are observed 
only for younger NS, and radiation from the young NS or 
the parent supernova pre-ionizes the ISM, suppressing Ha 
emission, or alternatively, that these younger NS mc more 
likely to be found in superbubbles which have insignificant 
neutral fractions. Future observations should resolve this 
issue. Meanwhile, it is satisfying to recognize that the 
scaling argument for the stand-off radius developed for Ha 
bow shocks can be applied to radio and X-ray detected 
synchrotron nebulae as well, in cases where a cometary 
morphology and a known radio pulsar indicate that a bow 
shock mechanism is at work. 

There arc several such radio pulsars which have asso- 
ciated plerionic bow shock nebulae at radio wavelengths, 
but not in Ha: for example, B1951-I-32 (Hester & Kulka- 
rni 1988), B1757-24 (the Duck, Frail & Kulkarni 1991), 
B1853+01 (Frail et al. 1996) and B1643-43 (Giacani et 
al. 2001). Each of these has a clear association with a 
supernova remnant, attesting to the relative youth of the 
object. The best available parameters for these pulsars 
and the stand-off angles of the associated bow shocks (ei- 
ther quoted or estimated from the radio images) are listed 
in Table 3. 

None of the listed pulsars has a parallax measurement, 
forcing us to rely on DM distances or estimates of the dis- 
tance to the associated supernova remnant for D. The 
velocities are estimated from a recent proper motion mea- 
surement for B1951+32 (Migliazzo et al. 2002) and from 
the NS age and displacement from the supernova rem- 
nant center for B1853+01 and B1643— 43. Ad hoc errors 
^ 25% are assigned to the latter estimates, but it is worth 
emphasizing that they may be larger. For example, for 
B1757— 24, a proper motion of 63-80 mas yr~^ has been 
estimated from its spindown age and displacement from 
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Table 3 





Radio 


A\D X-R,\Y Bow 


Shock Pad 


AMETERS 






NS 






V 

V 












(kpc) 


(km s-i) 


(erg s-i) 


{") 




B1643-43 


G341.2+G.9 


6.9 ± 2.1 


475 


35.55 


< 25 


1 


B1757-24 


G5.4-1.2 


4.6 ± 1.4 


< 550 


36.41 


1.5 ±0.4 


2,3 


B1853+01 


W 44 


3.3 ± 1.0 


- 375 


35.63 


< 1 


4 


B1951+32 


CTB 80 


2.4 ± 0.6 


280 ± 80 


36.57 


3.0 ± 0.8 


5,6,7 


B0906-49 




6.6 ±2.0 


~ 60 


35.69 


< 1 


8 


B1823-13 




4.1 ±1.0 




36.45 


< 40 


9 


B1929+10 




0.33 ±0.01 


163 ±5 


33.59 




10,11,12 


J0537-6910 


N157B 


~ 50 (LMC) 


~ 600 


38.68 


1.2 ±0.4 


13 


J0617.0+2221 


IC 443 


1.5 ±0.6 


- 250 


36 ±0.3 


8.5 ±2.0 


14,t 


The Mouse 


G359.2-0.5? 


5.0 ± 1.0 


~ 180 


36.3 ±0.3 


15 ±4 


15,16,t 



References.— (1) Giacani et al. (2001); (2) Frail & Kulkarni (1991): (3) ^ < 25 mas yr-\ 
Gaensler & Frail (2000); (4) Frail et al. (1996); (5) Strom (1987); (6) Migliazzo et al. (2002); 
(7) Hester & Kulkarni (1988); (8) Gaensler et al. (1998); (9) Finley, Srinivasan & Park (1996); 
(10) Yancopoulos, Hamilton & Helfand (1994); (11) Wang, Li & Begelman (1993); (12) Brisken 
(2001); (13) Wang et al. (2001); (14) Olbert et al. (2001); (15) Yusef-Zadeh & Bally (1987); 
(16) Predehl & Kulkarni (1995); (f: Not detected as a radio or X-ray pulsar.) 



the alleged center of the remnant G5.4— 1.2, but Gaensler 
& Frail (2000) find a 5(t upper limit < 25 mas yr^^ on its 
proper motion from interferometric observations. Only fu- 
ture interferometric measurements of the proper motions 
and parallaxes of these pulsars will reduce the large un- 
certainty associated with these estimates. 

As before, the estimated stand-off angles (or limits) arc 
plotted against a combination of observable parameters, 
and lines of constant density and inclination angle to the 
LOS are superimposed on the plot. It is immediately ap- 
parent that the large errors on the measurable parameters 
limit the useful constraints on the ISM density that can 
be extracted from the plot, and to illustrate the associ- 
ated uncertainties, the location of B1757— 24 using the er- 
roneous proper motion estimate of 63 80 mas yr"-'^ is also 
shown. However, it is reassuring that in spite of the large 
uncertainties, all the objects are consistent with physically 
plausible values of the ISM density. 

An unusual radio nebula has been observed for B0906— 49 
(Gaensler et al. 1998), which appears to be a pulsar wind 
nebula generated by a slow-moving (~ 60 km s~^) pulsar 
and confined by the ISM with a stand-off angle < 1". Us- 
ing an analysis similar to this work, Gaensler et al. (1998) 
conclude that a dense ISM is required (riA > 2 cm~^). 

Two other pulsars may have associated X-ray bow shock 
nebulae: B1823-13 and B1929±10. Both of these nebu- 
lae are ROSAT detections. For B1823-13, Finley, Srini- 
vasan & Park (1996) suggest that the observed extended 
component may be cither the associated supernova rem- 
nant for this young pulsar (spindown age ~ 21 kyr) or a 
bow shock nebula. If it is a bow shock nebula, the large 
scale size observed for the nebula (^ 40") implies either 
a low density medium or a low NS velocity. In the ab- 



sence of a proper motion estimate, the situation cannot 
bo resolved. For B1929+10, ROSAT data show a diffuse 
X-ray trail aligned with the pulsar proper motion, which 
has been explained as synchrotron radiation from a bow 
shock (Yancopoulos, Hamilton & Helfand 1994; Wang, Li 
& Begelman 1993). While an improved interferometric de- 
termination of the proper motion and parallax is available 
for this pulsar (Brisken 2001), the stand-off angle cannot 
be determined from the ROSAT data. Future observations 
will complete the picture for these two pulsars. 

In addition to synchrotron nebulae associated with 
known radio pulsars, there are two nebulae which are 
plausibly bow shocks associated with radio-quiet NS, 
the Mouse and the cometary nebula in IC 443, and 
a bow shock nebula associated with the X-ray pulsar 
J0537-6910. in the LMC. 

G359. 23-0.82 (the "Mouse") was discovered in VLA ob- 
servations by Yusef-Zadeh & Bally (1987). It appears to 
be a cometary wake, possibly emerging from the circular 
SNR G359.2— 0.5. Though no associated radio pulsar was 
known, Predehl & Kulkarni (1995) found an X-ray coun- 
terpart, and suggested that the object was most likely a 
young NS bow shock, analogous to the one produced by 
B1951+32. They also estimate «rE ^ Xif^'^ erg s^^ from 
the ROSAT X-ray spectrum. Yusef-Zadeh & Bally (1987) 
estimate a velocity of ~ 300 km s~^ for D « 8.5 kpc, which 
we scale down to 180 km s~^ for D w 5.0±1.0 kpc (Predehl 
& Kulkarni 1995). For a stand-off angle ~ 15", this im- 
plies an upper limit on the ISM density, n\ < 0.01 cm~^, 
which is fairly low. However, the previous caveats apply 
here, especially to the estimates of velocity and E. A faint 
radio pulsar counterpart has recently been discovered (F. 
Camilio, private communication), opening up the possibil- 



Bow Shocks from Neutron Stars 



11 




Fig. 7. — Scaling of stand-off angle with other observable param- 
eters for X-ray and radio-detected bow shock nebulae, as listed in 
Table 3. Circles mark radio pulsars and crosses mark the X-ray pul- 
sar J0537-6910, the radio-quiet NS CXOU J0617.0-I-2221, and the 
Mouse. The solid square shows B1757— 24 when its proper motion 
is (erroneously) estimated from its spindown age and displacement 
from the alleged supernova remnant center, emphasizing the possi- 
ble errors on this plot. Lines of constant number density n\ = 0.01, 
0.1, 1.0 cm~* are overplotted, for inclination to the line of sight 
i = 90° (thinnest to thickest solid lines) and i = 45° (thinnest to 
thickest dashed lines). 

ity that these estimates may be improved in future. 

A similar example is provided by the cometary nebula 
discovered by Olbert et al. (2001) in the supernova rem- 
nant IC 443. The NS producing this bow shock was de- 
tected by Chandra (CXOU J0617.0+2221), along with an 
X-ray counterpart to the radio synchrotron nebula. Based 
on the X-ray luminosity, E IQ^^-^ erg s~^. The distance 
to the associated SNR (and hence the NS) is ~ 1.5 kpc. 
From the bow shock morphology and the ratio of the NS 
velocity to the shock velocity in the SNR, Olbert et al. 
(2001) estimate a velocity of 250 ± 50 km s~^, indepen- 
dent of the distance. 

J0537— 6910 is a young Crab-hke X-ray pulsar (16 ms 
period) in the LMC, associated with the SNR N157B, 
which was serendipitously discovered in RXTE observa- 
tions (Marshall et al. 1998), but is not detected as a ra- 
dio pulsar. Chandra observations show large scale diffuse 
emission, which is again consistent with a ram-pressure 
confined pulsar wind nebula (Wang et al. 2001). The pul- 
sar has a measured spindown E = 4.8 x 10^^ erg s~^, and a 
velocity ^ 600 km s"-'^ estimated from its age and displace- 
ment from the SNR center. The estimated stand-off angle 
= 1'.'2 ± 0'.'4 is consistent with the bow shock description. 

The parameters for these objects are included in Table 3, 
and they arc plotted in Fig. 7. The uncertainties arc large, 
given the low confidence estimates oi D, E and velocity. 
However, constraints on the density of the ISM can still 
be obtained, and their location on the plot suggests that 



the bow shock model is a suitable one for these objects. 

6. DISCUSSION AND FUTURE WORK 

Bow shock nebulae, produced by NS moving supersoni- 
cally through the ISM, constitute unique probes of the en- 
ergetic environment where NS relativistic winds interact 
with the surrounding ISM. In this work, we report multi- 
wavelength observations of the Guitar Nebula, a spectacu- 
lar Ha bow shock nebula produced by a high velocity, mod- 
est E radio pulsar. High resolution HST observations show 
the detailed structure of the bright head of the nebula, 
confirming its limb-brightened nature and demonstrating 
the existence of a signific;ant asymnictric component. The 
complex shape is best explained by the existence of a den- 
sity gradient in the medium which is not aligned with the 
NS velocity vector. 

A simple momentum conserving bow shock model is 
used to fit the HST observations and extract the stand-off 
angle as well as placing joint constraints on the inclina- 
tion to the LOS, the density of the ambient medium, and 
the velocity and distance of the pulsar. While the low S/N 
ratio of the current HST image precludes a more detailed 
analysis, future HST data will provide more definitive con- 
straints for some of these parameters. 

The stand-off angle 9q and the Ha flux scale with 
other measurable parameters {E, D, jj) and unknowns 
(riH, X). These scaling relationships are derived from 
the basic physics of the bow shock formation. Using the 
observations of the Guitar Nebula (described here) along 
with published parameters for the other known Ha bow 
shocks, we confirm the scaling of 6*0 with the other pa- 
rameters for a diverse collection of nebulae, from simple 
canonical bow shocks to contorted shapes. The scaling 
relationship is used to derive upper limits on the density 
of the ambient medium, independent of the inclination of 
the bow shock to the LOS (which is hard to measure in 
a model-independent way). Additionally, the scaling of 
Ha flux is demonstrated with different combinations of 
measurable parameters: while the relationship is not as 
constraining in this case, the observations conform to the 
expected scaling within the inherent uncertainties in the 
neutral fraction and extinction and the unknown density 
of the ambient medium. 

We have also consolidated the available parameters for 
radio and X-ray bow shocks described in the literature, 
both for standard radio pulsars and radio-quiet or un- 
detected NS. The stand-off angle appears to scale as ex- 
pected for these objects as well. However, as described by 
Gaensler & Frail (2000) and illustrated by Fig. 7, the infor- 
mation available for these NS is inadequate, and has large 
margins of uncertainty. VLBI measurements of the paral- 
laxes and proper motions of these objects will be essential 
for a clearer understanding of the physics of plerionic bow 
shock nebulae and their correlation with young NS and 
well-deflned supernova remnants. 

The observed anti-correlation of Ha and plerionic radio 
emission from bow shocks is suggestive, though it needs to 
be confirmed in a larger sample of objects. Some young 
NS may be inside superbubbles created by previous su- 
pernovae, where the neutral fraction is insignificant. Pre- 
ionization of the medium, either by the young NS, or by 
the supernova in which it was born, may also suppress Ha 
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emission by reducing tlie cnutral fraction, leading to the 
prediction that young NS with associated supernova rem- 
nants should not produce Ha bowshocks. Conversely, if 
radiation from a young NS pre-ionizes the medium. Ha 
bow shocks should only be associated with older NS. This 
would be consistent with the apparent drop off in the ef- 
ficiency with which E is converted to radio emission in 
pulsar wind nebulae as a function of age (Gaensler et al. 
2000). 

Thus there arc multiple avenues of future enquiry that 
promise to be fruitful. As the recently reported discov- 
ery of a bow shock nebula for B0740— 28 demonstrates 
(Jones, Stappers & Gaensler 2002), deep observations of 
suitably selected pulsars will probably result in the detec- 
tion of more Ha bow shocks, leading to constraints on the 
NS and ISM properties, as described in this work. Fu- 
ture radio and X-ray observations will also produce more 
such objects, and clarify the relationship between the NS 
age and environment and the type of bow shock emission. 
Additionally, further interferometric proper motions (and 
parallaxes, where possible) are required to produce firmer 
constraints on the properties of the ISM and the NS using 
known bow shock nebulae, while theoretical models of bow 
shock evolution can elucidate the differences in the behav- 
ior of bow shocks downstream from the stand-off region. 
It is also possible that old high- velocity neutron stars, in- 
cluding some born in the Galactic halo, may be identified 
through the serendipitous discovery of radio, X-ray or Ha 
bow shock nebulae. 

The Guitar Nebula, specifically, will continue to be of 
interest, both from a theoretical modeling perspective, and 
for observations of its evolution over time. Monitoring of 
the pulsar dispersion measure and changes in the shape of 
the head of the Guitar in multi-epoch high resolution ob- 
servations will provide a unique perspective on the small- 
scale inhomogeneities in the interstellar medium. 
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